BCS-BEC Crossover in Mix-dimensional Fermi Gases 
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We investigate a mix-dimensional Fermi-Fermi mixture in which one species is confined in two- 
dimensional(2D) space while the other is free in three-dimensional space(3D). We determine the 
superfluid transition temperature T c for the entire BCS-BEC crossover including the important 
effects of noncondensed pairs. We find that the transition temperature reduces while the imbalance 
of mass is increased or lattice constant (dz) is reduced. In population imbalance case, the stability 
of superfluid is sharply destroyed by increasing the polarization. 
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1. INTRODUCTION 

Ultracold Fermi gases have attracted considerable attention [1-7] with their various tunability like the interaction 
strength via the Feshbach resonance and the dimensions of the system by means of optical lattices. This provides 
a wonderful opportunity to simulate and study the superfluidity of the fermions for entire BCS-BEC crossover. 
Lately, the population imbalance between the two species can also be varied. Since the imbalance causing the Fermi 
surface mismatch destabilizes pairing between the two species and enriches the ground state of the system [3, 4, 8- 
12, 14, 15]. The imbalance system possess exotic superfluid phases such as Sarma phase [14] mixture of the normal and 
superfluid phase, Larkin-Ovchinnikov-Fuldc-Fcrrcl (LOFF) state [12, 13] in which the condensed pairs have nonzero 
net momentum. 

Recently, the mix-dimensional system where the different species live in different spatial dimensions [16-19] has 
been realized with Bose-Bose mixture [20] by exploiting a species-selective ID optical lattice to confine only one 
atomic species( 41 i<') in 2D while the other ( 87 Rb) is free in the 3D space. In the mix-dimensional system, the Efimov 
effect [21-23] takes place only in some range of mass ratio for Fermi-Fermi mixtures [1G] while for any mass ratio for 
Bose-Bose and Bose-Fermi mixtures at resonance. The mix-dimensional systems open a new subfield in cold atom for 
investigating the heteronuclear molecules whose two species constituents live in different dimensions. 

In this paper we consider the two-species Fermi-Fermi mixtures in the mix-dimensional(2D-3D) system and obtain 
the behavior of the superfluid transition temperatures T c for the entire BCS-BEC crossover with different optical 
lattice parameters, the intra-specics mass ratio and polarization. In addition we give the behavior of T c in resonance 
limit. We consider only a uniformly superfluid of the mixtures excluding from consideration the phase separated state. 
As for the two-species living in different spatial dimensions, the Fermi surfaces of the two-species are mismatched 
which also can be tuned by the parameters of the optical lattice, mass ratio and the polarizations. We analyze the 
superfluid of the mix-dimensional system with mismatched Fermi surfaces including the the effects of noncondensed 
pairs which is called "pseudogap effects" [24, 25] based on the BCS-Leggett ground state. For the nonzero temperature, 
the excitation gap(A) contains two contributions from pseudogap for the noncondensed pairs(A pg ) and the order 
parameter(A sc ). The pseudogap contribution vanishes at zero temperature while the order parameter vanish at 
the critical transition temperature T c . The main conclusions are as follows: (a) increasing of the optical lattice 
parameter dz (half wave length of the optical lattice) can improve the transition temperature; (b) the masses and 
the popularizations imbalances of the two-species both suppress the transition temperature T c ; (c) in population 
imbalance case, the superfluid phase is in a narrow region(— 0.2 < p < 0.2) and unstable for high polarization. 

This paper is organized as follows. In Sec. II, the formalism of the mix-dimensional fermi gases are introduced. In 
Sec. Ill, the numerical results and discussions are given. The conclusion are given in Sec. IV. 
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2. FORMALISM OF THE MIX-DIMENSIONAL FERMI GASES 



The Hamiltonian for the mix-dimensional Fermi-Fermi mixtures system can be written as(S = ks = 1) 

k,<T k,k',q 

where, the pseudo-spin a = f , 4- labels the two types of Fermi atoms. We assume that the f Fermi atoms is confined 
in a species-selective ID optical lattice while the other is free in the 3D. The free dispersions of the two species Fermi 

k 2 -f/e 2 k 2 + fc 2 +/c 2 

atoms are £k,t = ; m " + 2i [1 — cos(fc z d z )] — /i^ and ^4 = " — 1 — Mi- Here is the two species Fermi atoms 
chemical potential, t is the amplitude of the t Fermi atoms tunnelling to the nearest-neighbor sites, and d z is the 
lattice parameter which is half wave length of the ID optical lattice. 

We determine the transition temperature T c including the pair fluctuation (pscudogap effect). Truncating the 
equation of motion for Green's functions scheme we can get the pair propagator 

1 + 9X(Q) 

and the self-energy 

E a (K) = Yl KQ)Go,a(Q ~ K \ (3) 
Q 

where the pair susceptibility is given by 

x(Q) = \lxn(Q) + Xit(Q)] = \ Y} G ^ ~ + G<u(Q - K)G t (K)}, (4) 

K 

with the bare green's function Gq a (K) = iu n — £k,er, and the G a (K) is dressed green's function. We take the notation 
K = (iw„,k), Q = (if2 n ,q), = rE n Ek, etc., where uj n (Q n ) is the odd(even) Matsubara frequency. 

Below the critical temperature T c , the T matrix and self energy contains both the condcnscd(sc) and the noncon- 
densed or "pseudogap"-associatcd(pg) contributions: 

t(Q)=t sc (Q) + t pg (Q), (5) 
t sc (Q) = -^S(Q), (6) 

UQ) = TT^(Q)' Q ^°- (7) 

So the total fermion self-energy is given by 

Z a (K) = ^sc(K) + ^a,pg(K) = -A 2 G . W (-K), (8) 

where the excitation gap contain two contributions of the order parameter (sc) and the pseudogap(pp) A 2 = A 2 c +A 2 g , 
with the pseudogap A 2 9 = - £q^o t pg (Q). 

The dressed green's function G a can be derived by Dyson equation 

!W " 4k «w„ + 

G ^)=^V+^4^ 5 do) 

l^n ~ Ck " + ^k 



where, the fermion excitation ^ = £ k ±4* k , with E k = + A 2 , £± = (£ k ,t±6u)A w 2 = (1 + £+/£; k )/2, « 2 

(l-£/23k)/2. 
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In supcrnuid state, the "gap equation" is given by the Thouless criterion, t _1 (Q = 0) = 0, which is equivalent to 
the BEC condition of pairs, /i pa ; r = 0. The coupling constant g is written in terms of the s-wave scattering length 
a, via the relationship M/(4na) = 1/g + £k(2e£), where M = 2m^m±/(m^ + mjj and e ± = (e-|- ± £f)/2. The gap 
equation reduce to 

-^s wl t fli! ' 4i (id 

k k £ k 

where the f(x) is the Fermi distribution function. The number of a fermion: n a = T,KG a (K). The equations of the 
total number n = n-j- + nj_ and the number difference pn — — of fcrmions are 

n = + + /(££))], (12) 



P" = £[/(£) -/(££)], (13) 

k 

here, p is polarization of the two-species Fermi atoms. Since the contribution to pseudogap is dominated by the small 
Q divergent region for the BEC condition. The T matrix can be expanded in low energy and long wavelength limit. 
The equation of pseudogap can be written as 

A L = -E^(Q) = ^ 1 EW' ( 14 ) 

Q q 

2 2 

where, fi q = + is pair dispersion with and M* representing the anisotropic effective pair mass computed 

from the low energy and long wavelength expansion of the pair susceptibility x(Q) [given in Appendix] and fo(f2 q ) is 
the boson distribution. 

The stability of the supcrnuid phase requires the positive defmiteness of the number susceptibility matrix, which 
is equivalent to the positive second-order partial derivation of the thermodynamical potential with respect to the 
excitation gap (c^O/dA 2 )^^ > 

, v a 2 1 -/(££)-/(££) /'(ff) + n 

k K 

In addition, we also consider the positivity of the anisotropic effective pair mass M* and Mj_ . Eqs.(ll)-(14) are closed 
for our mix-dimensional system. For T = 0, the pseudogap contribution vanish so the excitation gap (A = A sc ). For 
T > T c , the order parameter (A sc ) vanish such that the excitation gap contains the pseudogap only (A = A pg ). The 
transition temperature T c can be determined by sclf-consistcntly solving this closed equations throughout the entire 
BCS-BEC crossover with A 2 SC = 0. 



3. NUMERICAL RESULTS AND DISCUSSIONS 



First, the transition temperature T c , arc obtained from the self-consistent equations, as a function of scattering 
length, l/kpa, for different optical lattice parameter, kpdz, with the tunnelling amplitude t = Ep in population 
symmetry case (p = 0), shown in Fig.[l]. In the BCS limit, l/kpa — > — oo, the transition temperature T c as well as 
the anisotropic effective mass M z and M± of the cooper pairs is reduce to zero. As the interaction increase to the 
unitarity, T c rises rapidly. In the BEC limit, l/kpa — > oo, T c varies very slowly and the anisotropic effective mass 
of the pairs = 2m while M* depends on the value of kpdz. Increasing the parameter dz, lattice constant of the 
ID optical lattice, improves the transition temperature for the entire BCS-BEC crossover. This phenomenon can be 
illustrated as follows. When the tunnelling amplitude t = Ep, the Fermi surface of the J. is a spherical shell while the 
t is a prolate spheroid. Increasing of the parameter kpdz can decrease the mismatch of the two Fermi surface such 
that increasing kpdz is propitious to form the cooper pairs. If the parameter dz is much larger than the inverse of 
Fermi momentum kp 1 , the localization of the t fermion in z direction is very weak and the Mix-dimensional effect is 
reduced. Hence the parameter dz must not be much larger than the inverse of Fermi momentum kp 1 . 
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FIG. 1: (color online). T c as a function of 1/kpa for different lattice parameter kpdz with the polarization p — 0; t — Ef = 
h 2 k F /2m and = m;, Here fc_p is the noninteracting Fermi momentum for polarization p — 0. 




FIG. 2: (color online). Self-consistent solutions of the transition temperature T c as a function of l/fc^a for different mass ratio 
r\ = m-f/m_i with the polarization p = 0; t = -Ef = h 2 k%/2m and fcpdz = 1. The transition temperature T c reduce as the mass 
imbalance increasing. 



Second, the transition temperature T c , are obtained from the self-consistent equations, as a function of 1/fcpa for 
different mass ratio ry, shown in fig. [2]. We consider the population symmetry, the tunnelling amplitude t = Ep and 
the parameter kpdz = 1. In terms of the mixed dimensions system of Fermi-Fermi mixture(2D-3D), the system is 
stable against the Efimov effect [16, 22, 23] for even-parity case. 

Then, the transition temperature T c attains the maximum at the mass balance case and reduces for the mass ratio 
deviating from the symmetry case as shown in the Fig. [3] at resonance limit. The imbalance of the mass disturbs the 
inter-species pairing. 

Finally, for the population imbalance case, we find that the superfluid phase is restricted in the region of low 
polarization —0.2 < p < 0.2 at resonance, shown in Fig. [4]. The transition temperature T c reduces while the polar- 
ization is increased. The superfluid phase in the BEC limit is not stable for the condition of stability is violated as 
(d 2 il/dA 2 )^ < 0. At high polarization, the superfluid is unstable in the entire BCS-BEC region which is very 
different from the pure three dimensions systems [20]. The superfluid of the Fermi-Fermi mixture system in mixed 
dimensions is very unstable with population imbalance. In BCS side the anisotropic effective mass smoothly reduce 
to zero. In the unstable regions, other possible phase like the LOFF state maybe occur which will be consider in the 
future works. 
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FIG. 3: T c as a function of the mass ratio rj = m^/m\, at resonance 1/kpa = with the polarization p — 0; t = Ef = h 2 k F /2AI 
and kpdz — 1. 
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FIG. 4: T c as a function of the polarization p = (n-f — n±)/n at resonance 1/kpa — in the case of t = -Ef = h 2 kp/2m; 
kFdz = 1 and m-f = mj,. 

4. CONCLUSION 

To summarize, we obtain the transition temperature of the Fermi-Fermi mixture in Mix-dimensional system, as a 
function of the interaction strength for different optical lattice parameter dz and mass ratio in the BCS-BEC crossover. 
In addition, we also obtain the transition temperature behavior in resonance for different mass ratio and polarization. 
We obtain the transition temperature include the effect of the noncondensed cooper pairs by the generalized mean 
field theory based on the BCS-Leggett ground state in GqG schema. Increasing of the optical lattice parameter dz 
have positive effects to the transition temperature. In BEC limit, the anisotropic effective mass M± = 2m for different 
dz. In mass case, the transition temperature has maximum value at mass symmetry and reduce while the masses 
deviate from the balance. Superfluid phase is restricted within a narrow region of low polarization in population 
imbalance case. 
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Appendix: 

In this appendix, we evaluate the expansion of the pair susceptibility in low energy and long wavelength limit and 
derive the explicit expression of the anisotropic effective mass of the cooper pairs. 
The pair susceptibility is given by 

x(Q) = kxn(Q) + Xit(Q)] = \ E^tW - K)G X {K) + G 04 (Q - K)G f (K)], (16) 

Z Z K 

The explicit expression of the pair susceptibility can be got by substitute the Green's functions 

v(n) ly r l-/(ffi-/(&- q ,t) 2 /(ff)-/(&- q ,t) 2 
X V — o — 7a 7 7\ k — 7S — 7 1 — o - k 

2 k £k + &-q,t - if2 « ^k - 6t-q,t + ^ 



l-/(^)-/(fk-q4) 2 /(^)-/(fk-q,4), 



(17) 



At transition temperature T c , tJ^(Q) = and 

i PS 1 (Q) = i- 1 (Q) = x(Q)-x(o), (18) 

In low energy and long wavelength limit 

t-g(Q) = Z(iVL n - flq + fi pair + iT qt a), (19) 

2 2 

^q = 2m* + Tm* ^ s ^ ne dispersion of the cooper pairs. M z and Mj_ is the anisotropic effective mass. Below 
the critical temperature T c the chemical potential of cooper pair vanish fi pa i r = as for the BEC condition. The 
imaginary part of the pair dispersion is given by 



r q , Q = J - f(£) - /(6c- q ,tK*(££ + a- q .t - n) + [/(£) - /(&- q ,t)]^(« - &-q,t + n) 

k 

+ [1 - /(£) - /(ek-q,J]^*(fk-q4 - 0) + [/(d) - /(&-qJ]^(e£ - &-q4 + ( 2 0) 



and the reverse of the residue is 

dx 



2A 2 

The anisotropic effective masses of the cooper pair are given by 



din |n=q=0 

r n-S k (/(Ck,t) + /(ek4))], (21) 



1 1 d% 9 



2M* 4Z <9q 



2+-1 

— |f2=q=0 



1 v r 2 / , (^,t) / ^k,t > | 2 , %a 2 

8Z ^ A 2 1 dk ± j A 2 1 dk ± ' 

k 

(i - n€) mime* &) + (/(eg) - n&m + £) . g^t ) 

2£ k A 2 1 9ki J 

2£ k A 2 1 aki j 

(i - /(ffi - /(^))(^k - & 2 - urn - /(^))(g k + sir d^ 2 

£ k A 4 1 <9k ± ; 

(1 - /(g) - /(£))(£* - &) 2 - (/(d) - /(d))(^ + &) 2 , 36a , 21 

£ k A 4 1 <9k ± J 
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2M* 



2+-X 



1 d 2 t 
2Z~dkt~ |o=q=0 



k 

(i - fig) f(Q)(E k - C) + (/(eg) - /(e k ))(^ k + C) 9 2 6 



A 2 



1 dfc* ' 



A 2 



1 fljfc ' 



k,f 



2£ k A 2 9fc2 
(1 - /(Cg) - /(£))(13k - + - /(e k ))(^k + £) 



2E k A 2 



dk 2 



(i - f{£) ml))(E k - & 2 (/(g g)-/(^))(j?k+e k h ) 2 ( ^k, t , 2 



S k A 4 



9fe 



-)- 



(i - /(eg) - /(e k ))(£ k - e,[) 2 - (f( &- f(&)(E k +& 2 { d&, U2 



S k A4 



9^ 



) 2 }, 
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